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The  inv<  «tigations  reported  in  the  2nd  Semi-Annual  Technical 
Report  (Report  No.  9)  were  continued  during  the  above-mentioned  period 


in  the  following  directions: 

1.  strain  release  studies  of  the  circum-Pacific  seismic  belt) 

2.  laboratory  studies  of  the  stress  distribution  in  models  of 


geological  bodies. 

Below  we  shall  present  the  aim  and  the  results  obtained  in  the  reported 
time  interval  in  these  directions  of  our  investigations  of  the  rheologic 


properties  of  the  solid  earth. 


1.  Strain  release  studies 
a.  Aim 

The  aim  was  mainly  the  same  as  that  mentioned  in  the  2nd  Semi- 
Annual  Technical  Report  (Report  No.  9),  i.e.  to  increase  our  knowledge 
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of  the  dynamics  of  strain  release  during  earthquake  sequences,  and  the 
mutual  dependence  of  earthquake  activity  in  adjacent  areas.  We  included 
also  investigations  of  strain  release  patterns  preceding  large  earthquakes. 
The  emphasis  has  usually  been  concentrated  to  the  dynamics  of  aftershock 
sequences  in  earlier  studies.  However,  it  is  our  belief  that  much 
valuable  information  can  be  obteened  by  including  phenomena  preceding 
large  earthquakes.  In  addition,  a  world-wide  statistical  study  of  strain 
release  in  dependence  of  depth  was  carried  out  for  an  interval  of  35 
yearb,  the  aim  in  this  case  mainly  being  to  see  if  any  weak  layers  in  the 
mantle  would  exhibit  minima  in  the  strain-depth  curve. 

b.  Results 

The  paper  attached  to  the  2nd  Semi-Annual  Tech  ical  Report  with 
the  title*  Strain  release  in  the  circum-Pacific  belt*  Kern  County  1952, 

Des  rt  Hot  Springs  1946,  San  Francisoo  1957 »  by  Seweryn  J.  Duda  and 
Markus  B&th  has  been  accepted  for  publication  in  the  Geophysical  Journal 
of  the  Royal  Astronomical  Society,  London,  England,  and  is  at  present  in 
press. 

Two  additional  papers  have  been  prepared:  Strain  release  in  the 
circum-Pacific  belt*  Chile  I960,  by  Seweryn  J.  Duda,  and:  Strain  release 
in  relation  to  focal  depth,  by  Markus  BAth  and  Seweryn  J.  Duda.  The 
first  paper  i'  favored  for  publication  in  the  Journal  of  Geophysical 
Research,  American  Geophysical  Union,  Washington  D.C.,  USA;  the  second 
paper,  of  which  the  complete  manuscript  constitutes  an  Appendix  to  the 
present  report,  is  intended  for  publication  in  Geofisica  >ura  e  applicata, 
Milan,  Italy. 

The  study  of  strain  release  in  South  America  before  and  during  the 
Chilean  earthquake  sequence  starting  on  May  21,  I960,  supported  some 
older  results,  and  gave  some  new  ones.  A  distinct  concentration  of 
strain  release  was  found  towards  the  northern  and  southern  end  of  the 
aftershock  area.  The  segment  with  decreased  activity  in  between  was 
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slightly  displaced  in  the  interval  of  32  months  following  the  start  of 
the  sequence.  The  strain- release  characteristic  exhibits  two  branches 
with  different  shape,  clearly  separated  in  time.  The  strain  release  in 
the  aftershock  sequence  occurred  in  the  oscillatory  manner  previously 
found  for  a  number  of  aftershock  sequences. 

The  seismic  activity  in  South  America  during  40  months  preceding 
the  Chilean  earthquake  soquence  exhibited  a  repeated  northward  migration 
of  the  maximum  strain  release,  the  migration  velocity  increasing  with 
time  towards  the  start  of  the  sequence.  Pour  migration  cycles  could  be 
observed,  the  pattern  of  which  was  practically  destroyed  by  the  beginning 
of  the  Chilean  earthquake  sequence.  In  particular, the  region  adjacent  to 
the  aftershock  area  in  the  north  became  quiet  for  more  than  two  years 
after  the  start  of  the  sequence* 

In  the  Chilean  sequence  the  totally  accumulated  strain  was  released 
by  several  strong  (main)  shocks  and  a  relatively  unimportant  aftershock 
sequence,  from  the  energetic  point  of  view.  This  is  in  contradiction 
to  some  other  circum-Pacific  aftershock  sequences,  e.g.  Kamchatka  1952 
and  the  Aleutians  1957,  where  the  totally  accumulated  strain  was  released 
by  the  main  shock  and  a  powerful  aftershock  sequence. 

A  similar  investigation  of  strain  release  has  been  started  for 
Western  North  America  for  the  time  interval  1955-1959*  Tic  strain 
release  pattern  prior  to  the  strong  Alaskan  earthquake  of  July  10,  1958 1 
has  several  features  in  common  with  the  strain  release  in  Chilo  up  to 
May,  i960.  Thus,  three  migration  cycles  could  be  stated  in  the 
investigated  time  interval  before  the  Alaskan  earthquake.  The  migration 
velooity  increased  towards  the  start  of  the  sequence.  The  migration  is 
directed  away  from  the  region  with  highest  stress  concentration,  where 
later  the  important  strain  release  occurred.  The  initial  rupture  in  the 
main  shocks  was  both  in  Alaska  and  in  Chile  directed  oppositely  to  the 
migration  direction.  By  the  strong  Alaskan  shock  the  migration  pattern 
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was  also  destroyed.  A  quietness  could  be  observed  in  the  region 
adjacent  to  the  aftershock  area. 

This  study  is  not  yet  finished.  A  detailed  report  is  planned  to 
be  given  in  the  near  future. 

The  study  of  world-wide  strain  release  in  dependence  of  depth  for 
35  years  was  performed  with  the  purpose  to  see,  if  the  asthenosphere 
low-velocity  layer,  known  from  seismic  wave  propagation,  gives  rise  to 
a  smaller  amount  of  strain  release,  than  expected  from  a  regular  strain 
decay  with  depth  in  the  corresponding  depth  range.  The  strain  release 
was  found  to  be  highest  in  the  uppermost  75  km  of  the  earth.  It  decreases 
exponentially  with  depth  between  75  and  400  km.  After  a  pronounced  minimum 
between  400  and  475  km  it  increases  again  approximately  exponentially 
between  475  and  650  km.  Thereafter  it  drops  rapidly  to  ^ero.  There  is 
no  pronounced  minimum  corresponding  to  the  asthenosphere  low-velocity 
layer. 

The  depth  curve  for  the  number  cf  shocks  is  nearly  parallel  to 
the  strain-depth  curve.  The  average  strain  per  earthquake  shows  only 
an  insignificant  decrease  with  depth. 

As  continuation  of  this  part  of  our  work  it  is  planned  to  carry 
out  a  strain-release  investigation  for  one  further  circum-Pacific 
earthquake  sequence  -  the  Kurile  sequence,  which  started  on  November  6, 
1958»  and  afterwards  for  two  Asiatic  earthquake  sequences  -  the  Mongolian 
sequence  which  started  on  December  4»  1957 »  and  the  Chinghai  Province, 
China,  sequence  which  started  on  May  21,  1962.  Perhaps  some  additional 
smaller  sequences  will  be  included.  One  aspect  of  this  research  is  to 
compare  the  strain  released  in  the  main  shock  to  that  of  the  aftershock 
sequence.  We  have  already  seen  that  these  strains  are  different  in 
different  regions.  This  may  give  a  clue  to  some  of  the  dynamic 
characteristics  of  different  seismic  regions. 
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In  all  this  research,  Benioff’s  methods  have  been  of  basic  importance* 
However,  it  would  be  desirable  t-»  improve  the  methods  as  far  as  possible. 
Especially,  we  aim  at  a  relation  between  magnitude  of  a  shv.  k  and  the 
volume  of  strained  rock,  instead  of  assuming  this  volume  as  constant,  as 
hitherto  done  in  all  such  studies . 

2.  Laboratory  investigations 

a.  Aim 

Stress  distributions  in  models  of  geological  bodies  can  be  obtained 
in  an  analytical  form  only  in  the  case  of  simple  forms,  and  usually  under 
the  assumption,  that  the  material  deformed  by  the  stress  is  perfectly 
elastic*  On  the  other  hand,  earthquake  sequences,  as  well  as 
tectonophysical  processes  in  general,  cannot  be  considered  as  perfectly 
elastic  processes i  it  is  necessary  to  introduce  the  time  parameter  into 
the  stress-strain  relation.  This  gives  rise  to  essential  mathematical 
difficulties.  A  different  way  to  get  quantitative  information  about  the 
stress  field  and  its  changes  with  time  is  an  analogue  processing.  We 
decided  to  apply  the  photoelastic  method,  first  to  the  steady  stress 
state  in  perfect  elastic  bodies,  and  thereafter  in  bodies  with  creep 
features . 

b.  Results 

Our  collaboration  with  the  photoelastic  laboratory  of  the  Royal 
Water  Power  Board  (Kungl,  Vattenfallsstyrelsen),  Stockholm,. gave  the 
first  result.  A  model  of  a  crack  in  a  plate,  the  plate  being 
sufficiently  large  to  be  considered  as  an  infinite  plane,  was  prepared. 

It  was  subjected  to  a  uniaxial  compression  under  45°  to  the  crack  axis. 

The  isochromatic  fringes  were  obtained  in  the  polariscope  and  the 
principal  stresses  will  be  calculated  for  this  case.  This  experiment 
permits  determination  of  the  stress  distribution  in  the  presence  of 
shear  stress  at  the  crack.  On  the  other  hand,  the  analytical  expressions 
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for  the  stress  distribution,  as  given  by  C.E.  Inglis  (Stresses  in  a  plate 
due  to  the  presence  of  cracks  and  sharp  corners,  Trans.  Inst.  Naval 
Architects,  Vol.  55,  Pt.  1,  1915,  pp.  219-250),  and  K.  Wolf  (Beitrage 
zur  ebenen  Elastizit&tstheorie ,  Zeitschrift  f.  techn.  Physik,  Vol.  2, 

Mo.  8,  1921,  pp.  ,'09-2l6)  are  calculated  under  the  assumption,  that  tho 
crack  boundaries  are  free.  Our  picture  of  isochromatic  fringes  was 
attached  to  the  Ninth  Monthly  Status  Report  (Report  No.  12).  The 
experiments  can  easily  be  varied  from  pure  tangential  to  pure  normal 
stress  at  the  crack. 

The  experiments,  thus  far  performed,  give  the  stress  distribution 
in  the  static  case.  It  is  naturally  of  interest  to  enlarge  these 
experiments,  so  as  to  be  able  to  study  also  the  dynamic  processes  at  a 
crack,  resembling  the  earthquake  phenomena  as  much  as  possible,  including 
the  time  factor.  Different  models,  e.g.  with  zones  with  decreased 
strength,  as  well  as  varying  stress  fields  are  envisaged. 


The  third  item  of  our  investigations  of  the  rheologic  properties 
of  the  solid  earth,  i.e.  the  development  of  methods  to  measure  stress 
and  strain  variations  in  the  earth's  crust,  was  not  advanced  much  in  the 
reported  time.  Wo  have  again  consulted  Professor  Hast  in  Stockholm 
about  his  stress  measurements  and  hope  to  embark  on  these  problems  in 
more  detail  when  items  1  and  2  have  been  more  advanced  or  fin.' shed. 

__  _ 

Markus  BAth 
Project  Scientist 
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STRAIN  RELEASE  IN  RELATION  TO  FOCAL  DEPTH 
by  MARKUS  BATH  (X)  &  SBVfERYH  J.  DU  DA  (XX) 

Summary  -  The  world-wide  strain  release  in  relation  to  focal  depth 
has  been  calculated  for  all  shocks  with  magnitude  7  and  over  for  the 
interval  1918-1952.  The  strain  exhibits  a  strong  maximum  in  the  upper¬ 
most  75  km  of  the  earth;  it  decreases  exponentially  with  depth  between 
75  and  400  km,  with  an  unimportant  minimum  corresponding  co  the 
asthenosphere  low-velocity  layer  anc  another  minimum  at  275  km;  after  a 
pronounced  minimum  between  400  and  475  km  it  increases  again  approximately 
exponentially  between  475  and  65O  km,  after  whicn  it  d^ops  rapidly  to 
zero.  The  shape  of  the  strain-depth  curve  is  interpreted  in  terms  of 
the  physical  conditions  and  the  intensity  of  strain  accumulation.  In 
particular,  the  increase  between  475  and  65O  km  is  ascribed  to  a 
combined  effect  of  temperature  and  pressure  variation  with  depth  with 
related  phase  changes  and  possible  changes  in  composition.,  The  depth 
curve  for  the  number  of  shocks  is  nearly  parallel  to  the  strain-depth 
curve,  and  the  average  strain  per  earthquake  shows  only  an  insignificant 
decrease  with  depth. 

Zusammenfassung  -  Es  wird  die  Tiefenabhangigkeit  der  Deformations- 
ausloeung  in  alien  Erdbeben  mit  Magnitude  7  vnd  dariiber,  im  Zeitraum 
1918-1952  untersucht.  Die  Deforraationsauslosung  hat  ein  ausge3prcchenes 
Maximum  in  den  obersten  75  km.  Sie  nimmt  im  Tiefenbereich  zwischen  75 
und  400  km  exponential  mit  der  Tiefe  ab,  wobe?  sich  ein  Minimum,,  das  der 

(X)  Dr.  MARKUS  BATH,  Seismological  Institute,  Uppsala.  Sweden. 

(XX)  Mr.  3EUERYN  J.  DUDA,  Seismological  Institute,  Uppsala.  Sv/eden  . 

On  leave  from  3EI3M0S  G.m.b.H.,  Hanover.  Germany, 


Schicht  rait  nieariger  Wellengeschwindigkeit  entspricht,  schwach  andeutet, 
und  ein  zweites  Minimum  bei  275  km  liegt.  Nach  einem  ausgepragten 
Minimum  zwischen  400  und  475  km  Tiefe,  steigt  die  Deformationsausldsung 
erneut,  etwa  exponential,  im  Tiefenbereich  zwischen  475  und  650  km  an, 
und  nimmt  danach  schnell  ab.  Per  Verlauf  der  Deformationsausllsung  in 
Abhangigkeit' vor  der  Tiofe  wird  an  Hand  der  physikalischen  Verhaltnisse 
im  Frdkorper,  sovie  der  Inten3itat  der  Deforraationsaufspeicherung 
gedeutet.  Insbesondere  wird  der  Anstieg  zwischen  475  und  65O  km  einer 
Anderung  der  Temperatur  und  des  Druckes  mit  der  Tiefe,  bei  Phasensprung 
und  moglicher  Anderung  in  der  Zuoammensetzung  des  Materials  zugeschrieben. 
Die  Tiefenabhiingigkeit  der  Anzahl  der  Beben  und  der  Deformationsauslosung 
verlaufen  annahemd  parallel)  die  mittlere  Doformationsauslbsung  pro 
Erdbeben  zeigt  lediglich  eine  unbedeutende  Abnahme  mit  der  Tiefe. 

1.  Introduction  -  There  is  no  doubt  that  £ei3mological 
measurements  offer  the  most  reliable  information  about  the  physical 
conditions  of  the  earth's  deep  interior.  Determinations  of  travel  times 
of  seismic  waves  and  of  seismic  wave  velocities  in  the  earth's  interior 
constitute  the  most  important  contribution  df  seismology  in  this  respect. 
But  it  is  clear  that  amplitude  measurements  will  furnish  additional 
information  of  great  value.  For  instance,  atte;  uation  of  seismic  waves 
deperis  on  the  material  properties  and  it  is  often  related  to  the 
frequency  of  the  waves,  causing  a  selective  attenuation  with  ensuing 
changes  of  the  wave  spectrum  during  propagation. 

Another  applicationof  amplitude  measurements  is  made  in  determinations 
of  magnitudes  and  strain  release  in  earthquakes.  The  present  paper  ia 
concerned  with  this  problem,  in  particular  the  variation  of  the  strain 
release  with  depth  and  its  possible  relations  to  the  physical  conditions. 

A  starting-point  for  this  piece  of  research  was  the  question  if  any  weak 
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layers  with  increased  plastic  flow  would  show  up  as  minima  in  the  curve 
relating  strain  release  +o  focal  depth. 

Strain  characteristics  of  the  earth's  interior  have  been  studied 
earlier  by  several  other  au  hors  (see  e.g.  BENIOPP  &  GUTENBERG,  1951 » 
BENIOFF,  1954|  GUTENBERG,  1958  a,  1959).  Earthquake  energy  in  relation 
to  focal  depth  was  studied  by  GUTENBERG  (195^,  1957)*  To  our  knowledge, 
the  only  earlier  investigation  of  strain-depth  relations  has  been 
published  by  RITSEMA  (1954)*  The  present  paper,  however,  differs  from 
that  one  in  several  respects,  as  in  methods  used  and  in  the  discussion 
of  the  results. 


2.  Observed  depth  relation  of  strain  release  and  of  number  of 
earthquakes  -  The  material  listed  by  GUTENBERG  &  RICHTER  (1954)  was 
used.  In  order  to  get  a  homogeneous  set  of  data,  we  restricted 
ourselve3  to  earthquakes  with  magnitue  M  -  7  in  the  interval  1918-1952. 
No  separation  according  to  earthquake  regions  has  been  made,  as  we  only 
aimed  at  the  world  distribution  of  strain  release  with  depth. 

In  calculating  the  strain  release  S,  we  used  for  comparison 
purposes  the  same  formulae  as  in  earlier  similar  studies,  i.e.,  with  J  - 
seismic  wave  energy, 


extending  the  sum  over  all  earthquakes  of  given  depth  in  the  range  of 
time  and  magnitude  a-  already  defined,  and 


log  J®  =  4*5  +  0.9 


Only  logarithms  to  the  base  10  are  used  in  this  paper* 

Equation  (l)  implies  that  the  following  quantities  are  constant: 
(a)  the  fraction  of  potential  strain  energy  converted  into  seismic  wa/o 
energy?  (b)  the  elastic  parameters}  (e)  the  volume  of  strained  rock 
for  every  earthquake (compare  BATH  &  BENIOFF,  1958)*  The  elastic 
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parameters  certainly  vary  with  depth,  but  within  the  depth  range  of 
33-600  km  the  square  root  of  the  elastic  parameters  concerned  varies 
roughly  between  0.8  and  1.1.  For  this  reason,  assumption  of  constant 
elastic  parameters  is  very  well  permitted  in  a  first  approximation  and 
dees  not  at  all  change  the  more  important  trends  of  the  strain-depth 
curve.  With  regara  to  (c)  it  is  known  that  the  volume  increases  with 
the  magnitude,  but  no  very  reliable  relation  has  so  far  been  given,  and 

the  assumption  of  constant  volume  is  equally  permitted  in  a  first 

approximation.  Effects  of  the  factor  (a)  are  discussed  below. 

Fig.  1  shov-r  the  results  of  the  computations.  The  strain  as  well 
as  the  number  of  shocks  with  Mi1/  has  been  plotted  for  every  25  km  in 
depth,  i.e.  corresponding  to  the  approximate  accuracy  of  focal  depth 
determinations . 

The  strain-depth  curve  in  Fig.  1  exhibits  a  series  of  features 
for  different  ranges  of  the  depth  hs 

(a)  h  <  75  km*  The  depths  of  shocks  in  this  range  are  not  well  known 

and  given  only  sparsely  by  GUTENBERG  &  RICHTER  (1954/*  As  a  consequence 

no  reliable  curve  can  be  constructed  for  this  range.  The  total  strain 
in  this  range  is  known  with  the  same  accuracy  as  for  greater  depths, 
but  the  depth  variation  is  uncertain.  In  Fig.  1  a  logarithmic  decrease 
with  depth  is  tentatively  traced.  Anyway,  it  is  clear  that  the  strain 
release  far  exceeds  that  -f  any  greater  depths  and  decreases  with  depth 
more  rapidly  than  at  greater  depths.  It  is  likely  that  the  strain  has 
a  maximum  at  some  shallow  depth,  perhaps  around  25-30  km,  and  decreases 
to  a  low  value  closer  to  the  surface  of  the  earth  (h  *  0).  This 
interprets  on  would  agree  with  MATSUSHIMA'S  (1961)  results.  KONING's 
(1953)  frequency-depth  curve  exhibits  a  relative  minimum  at  60-80  km 
and  a  relative  maximum  at  120  km  depth,  which  is  probably  only  an 
effect  of  insufficient  knowledge  of  focal  depth,  especially  in  the 
shallower  rang^. 
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(b)  75  km  £  h  £  400  km.  The  logarithr  of  the  strain  shows  an 
approximately  linear  decrease  with  depth.  A  least-square  solution, 
illustrated  by  the  straight  line  X  in  Pig.  1,  yielded  the  following 
equation,  disregarding  the  proportionality  f-ccor  in  eq.  (l)t 

(3)  log  S  =  (2.85  ±  0.36^  -  (0.0040  +  O.OOIO)  h 

iO  4 

the  unit  for  S  being  10  (ergs)®.  There  are  3ome  superimposed 
variations  in  the  lower  part  of  this  range,  especially  a  minimum  at  h  = 
275  km.  However,  the  number  of  earthquakes  is  also  quite  small  in  this 
range,  and  greater  fluctuations  are  to  be  expected.  Local  influences 
are  then  greater,  and  this  minimum  may  for  instance  largely  depend  on 
the  depth  distribution  in  the  Japanese  area  (see  GREENBERG  &  RICHTER, 
1954). 

(c)  400  km  <  h  475  km.  There  is  no  earthquake  with  M  i  7  in  this 
range  and  consequently  no  eoi re upending  strain. 

(d)  475  km  ^  h  ^  650  km.  The  strain  exhibits  e  pronounced  increase 
with  depth,  approximated  here  by  the  straight  line  B  in  Fig.  1,  which 
corresponds  to  our  least-square  solution: 

(4)  log  S  =  (-1.68  +  0.28)  +  (0.0054  +  0.0017)  h 

(e)  h  >  650  km.  There  is  no  shock  with  M  2^7  in  this  range.  As  is  well 
known,  no  shock  deeper  than  720  km  has  been  r — orded,  and  the  strain 
disappears  completely  at  this  depth  and  below. 

Pig  1  shows  also  the  depth  variation  of  the  number  N  of  earthquakes 
with  H  >7»  This  curve  exhibits  almost  perfect  parallelism  with  the 
strain-depth  curve.  This  will  be  approximately  true  even  if  we  had 
chosen  a  lower  magnitude  limit,  because  of  the  statistics  given  by 
GUTENBERG  &  RICHTER  ( 1954 *  P»  17»  eqs  5  and  6)  for  the  earth  as  a  whole. 
This  means  that  the  variation  of  strain  with  depth  is  ma'nly  a 
consequence  of  the  variation  of  the  number  of  shocks  with  depth.  As  a 
result,  the  average  strain  release  per  earthquake,  i.e.  S/N,  also  given 
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in  Pig.  1,  shows  only  an  insignificant  decrease  with  depth.  The  least- 
square  solution  for  the  depth  range  40  km  ^  h  ^  650  km  is  shown  by  the 
straight  line  C  in  Fig.  1  and  has  the  <.  ration: 

(5)  S/N  =  (11.97  +  2.94)  -  (0.0037  +  0.0033)  h 
There  is  only  a  minor  decrease  of  the  maximum  magnitude  with  focal  depth, 
as  is  evident  from  the  following  table* 

h  km  100  150  200  250  300  350  400  500  550  600  650 

M  max  8.0  8.0  7.8  7.3  7.2  7.8  7.8  7.2  7*3  7.8  7.6 


3-  Physical  interpretation  of  the  strain-depth  curve  -  Even  if 
our  curves  probably  represent  fair  approximations  to  the  general  conditions 
within  the  earth,  it  must  be  bom  in  mind  that  they  are  strictly  true 
only  for  the  limited  interval  of  35  years  considered.  In  evaluating 
the  physical  implications  of  our  results,  we  have  to  conside.r  that  only 
strain  which  has  produced  elastic  waves,  is  involved  in  our  observations. 
This  strain  depends  on  two  main  factors: 

(a)  The  strain  accumulation.  Without  accumulation,  there  will  naturally 
be  no  release.  The  strain  accumulation  depends  on  3low  relative  motions 
in  the  earth's  interior,  the  reasons  for  which  we  do  not  need  to  consider 
in  this  connection. 

(b)  The  mode  of  strain  release.  Depending  on  their  physical  condition, 
the  rocks  are  able  to  store  larger  or  smaller  amounts  of  strain.  With 
increasing  depth,  plastic  flow  processes  will  become  more  and  more 
important.  Ultimately,  they  are  prevailing  and  no  (shearing)  strain 
can  be  stored  for  any  length  of  time  of  importance  for  earthquake 
generation. 

Obviously,  the  variation  of  strain  with  depth  depends  on  both 
factors,  (a)  and  (b),  which  both  vary  with  depth.  This  presents  a 
major  difficulty  in  interpreting  the  empirical  stre  .i-depth  curve.  If 


(a)  varies  but  (b)  is  constant,  the  frequency  of  shocks  N  will  vary  in 
parallelism  with  the  intensity  of  strain  accumulation,  whereas  the 
maximum  magnitude  will  be  essentially  constant.  Variations  of  magnitude 
require  variations  of  the  factor  (b).  Our  empirical  results  would  thus 
indicate  that  variations  of  (a)  rather  than  of  (b)  are  of  dominating 
influence  on  our  strain-depth  curve. 

The  case  with  constant  (a)  but  variable  (b)  can  be  better  grasped 
with  the  following  simple  consideration.  Compare  two  cases,  1  and  2, 
both  having  the  same  rate  of  strain  accumulation,  i.e.  the  same  (a), 
but  with  more  of  plastic  flow  processes  in  2  than  in  1,  i.e.  (b)  is 
different.  Denoting  the  two  cases  by  indices  1  and  2  and  applying  the 
equation  for  J  expressed  in  M,  we  have 

(6)  log  (NjJ^Jg)  -  t(M1  -  M2)  -  log  (N2/N1) 

with  *he  constant  b  =  1.8  according  to  eq.  (2)  but  •  1.4-1 .5  according 
to  more  recent  determinations.  In  eq.  (6)  we  have  for  simplicity 
assumed  constant  magnitudes:  in  case  1  and  2^  in  case  2,  corresponding 

to  the  breaking  strength  in  the  two  cases.  We  now  have  that 
because  of  the  lower  breaking  strength  in  case  2.  Moreover,  the  total 
strain  energy  is  the  same  in  the  two  cases,  as  (a)  is  the  same.  But, 
part  of  the  total  strain  is  dissipated  by  plastic  flow  in  case  2,  and 
therefore  the  fi action  of  total  strain  converted  into  elastic  wave 
energy  is  then  less,  i.e.  N^J^/NgJg  >  1«  Eq.  (6)  then  reduces  to  the 
following  condition: 

(7)  b(Mx  -  M2)  >  log  (N^) 

Thus,  may  be  greater  than  N^,  up  to  a  limit  defined  by  eq.  (7).  But 
it  is  not  necessary  that  exceeds  contrary  to  some  statements 
in  the  literature  that  a  zone  with  lower  breaking  strength  should  be 
accompa  .  -ed  with  a  larger  number  of  shocks. 

In  the  real  earth  the  cases  1  and  2  do  not  exist  isolated  but 


-  8 


may  be  adjacent  to  each  other,  which  will  modify  the  situation  to  some 
extent.  If  a  layering  corresponding  to  1-2-1,  i.e.  a  weak  layer  between 
two  stronger  layers,  is  subjected  to  a  shearing  strain,  it  will  first 
break  in  the  layer  2,  i.e.  the  whole  structure  will  roughly  behave  as 
only  composed  of  its  weakest  material.  This  means  that  the  effects  of  a 

*  i 

weak  layer  ext  and  upward  and  downward  probably  several  times  its  own 
thickness.  Beyond  a  certain  distance  the  weak  layer  has  no  longer  any 
effeot  on  the  strain  release. 

Penetrating  into  more  detail,  however,  the  conditions  are  rather 
complicated.  The  factors  (a)  and  (b)  cannot  be  considered  as  independent 
of  eaoh  other.  The  time  factor  is  of  importance,  and  in  fact  (b)  has  to 
be  taken  as  dependent  on  the  faotor  (a),  A  given  material  under  given 
physical  conditions  will  behave  differently  according  to  the  rate  of 
strain  accumulation. 

The  large  strain  release  within  the  crust  and  upper  mantle  down 
to  about  75  km  depth  is  explained  by  large  strain  accumhlation  (in  active 
earthquake  areas)  combined  with  large  ability  of  strain  storage. 

According  to  GUTENBERG  &  RICHTER  ( 1954) *  "isostasy  and  post-glacial 
uplift  indicate  that  the  strength  (resistance  to  plastic  flow)  below 
80+  kilometers  is  less  than  l/lOO  of  that  near  the  surface". 

The  existence  of  the  asthenosphere  low-velocity  layer  at  depths 
around  100-200  km  is  now  generally  recognized  from  wave  propagation 
studies.  It  is  considered  as  a  weak  zone,  where  one  shouxd  expect 
plastic  flow  processes  to  be  more  important.  The  zone  is  characterized 
by  minima  of  elastic  wav<-  velocities  and  of  elastic  parameters  (BATH, 
1956),  However,  our  strain-depth  curve  does  not  exhibit  any  clear 
minimum  to  be  ascribed  to  the  asthenosphere  low-velocity  liyer.  There 
is  only  a  small  superimposed  minimum  between  the  depths  100  km  and  225  km. 
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'  probable  reason  is  that  the  material  in  the  asthenosphere  still  has 
sufficient  strength  to  accumulate  considerable  amount  of  3train. 

Various  geophysical  effects  of  a  soft  or  weak  layer  in  the  upper  mantle 
have  been  discussed  by  SCHEIDEGGER  (i960),  HALES  U96l),  SHIMOZURU  (1963) 
and  others.  In  evaluating  effects  of  such  a  layer  it  is  naturally 
essential  to  consider  the  time  factor. 

The  strain  minimum  at  275  km  depth  could  possibly  be  due  to  a 
weak  layer  in  this  depth  range. 

The  most  significant  feature  of  the  strain-depth  curve  is  the 
pronounced  minimum  in  the  depth  range  of  400-475  km.  Our  curve  is  summed 
up  for  the  whole  earth.  In  many  earthquake  regions  there  are  no  deep 
shocks,  i.e.  the  secondary  maximum  below  475  km  does  not  occur.  The 
reason  is  probably  lack  of  strain  accumulation  at  these  depths  in  many 
regions.  The  accumulation  is  unevenly  distributed  over  the  surface  of 
the  earth,  and  its  variation  with  depth  is  most  likely  also  different  in 
different  areas.  On  the  other  hand,  the  physical  properties  at  these 
depthr  in  the  earth  are  probably  the  same  everywhere.  In  other  words, 
the  lack  of  deep  earthquakes  in  some  regions  is  rather  to  be  ascribed  to 
lack  of  strain  accumulation  than  to  plastic  flow  processes. 

It  remains  to  be  explained  why  the  material  behaves  in  such  a 
way  as  to  give  a  pronounced  strain  minimum  at  400-475  km  depth  and  a 
strong  increase  between  475  km  an:‘.  650  km.  BULLEN’s  (1947)  results 
indicate  a  second  order  discontinuity  at  a  depth  of  41?  km  with  rapid 
downward  increases  of  wave  velocities,  density  and  elastic  parameters. 
This  would  imply  increased  strength  of  the  material,  possibly  due  to 
phase  changes  or  changes  in  composition  (BIRCH,  I96I}  RINGW00D,  1962). 

For  the  interval  475-650  km  this  increase  of  strength  will  counteract 
the  general  tendency  towards  plastic  flow  with  depth  and  will  correspond 
to  a  greater  ability  to  store  strain.  Below  65O  km,  plastic  flow 
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processes  again  dominate  the  picture. 

The  strain  increase  between  475  km  and  650  km  may  thus  be  traced 
back  to  the  combined  effect  of  pressure,  temperature,  phase  changes 
and/or  changes  in  composition.  Increased  pressure  generally  means 
increased  strength,  whereas  increased  temperature  means  decreased  strength 
(MATSUSHIMA,  196l).  If  within  the  interval  mentioned,  the  pressure 
effect  dominates  over  the  temperature  effect,  this  would  explain  our 
results.  As  a  matter  of  fact,  the  ratio  of  pressure  increase  (dp)  to 
temperature  increase  (dT)  with  depth,  increases  rapidly  below  400  km 
depth.  Combining  BULLEN's  (1947)  pressure  curve  and  GUTENBERG'S  (1951) 
temperature  curve,  we  find  the  values  given  in  Table  1.  With  P  = 
breaking  strength,  the  condition  for  increasing  strength  without  change 
of  phase  or  composition  is* 

(8)  dp/dT  >  -  (*p/»T)p  t  (3P/?p)  T 
Unfortunately,  lack  of  knowledge  of  numerical  values  of  the  partial 
derivatives  under  the  extreme  conditions  at  great  depth  prevents  further 
calculations.  In  addition,  pressure  and  temperature  variations  with 
depth  will  determine  any  phase  changes  that  may  occur.  I£  the  phase 
change  possibly  in  combination  with  a  change  in  composition  is  spread  out 
over  a  depth  range  of  several  hundred  kilometers,  it  is  not  necessary  to 
have  a  second  order  discontinuity,  i.e.  the  velocity  curves  could  be 
more  in  agreement  with  those  of  GUTENBERG  (1958  b). 

On  the  other  hand,  as  emphasized  by  GUTENBERG  (see  BENIOFF  & 
GUTENBERG,  1951 ),  "aeep-focus  earthquakes  do  not  necessarily  require 
appreciable  strength  for  an  accumulation  of  stresses,  but  a  high 
coefficient  of  viscosity  is  sufficient  to  permit  the  required 
accumulation  of  strain". 

Explanations  in  terms  of  special  strain  accumulation  patterns  could 
also  be  imagined,  but  these  could  hardly  be  pushed  above  the  purely 
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hypothetical  level  (sea  e.g.  HESS,  1951). 

4*  Conclusions  -  The  world-wide  strain  release  in  relation  to 
focal  depth  for  all  earthquakes  with  M  >  7  in  the  interval  1916-1952 
has  revealed  the  following  features* 

(a)  The  maximum  strain  occurs  within  the  uppermost  75  km  of  the  earth, 
due  to  a  combination  of  strong  strain  accumulation  (in  seismically 
active  areas)  and  good  possibilities  of  strain  storage. 

(b)  The  asthenosphere  low-velocity  layer  pro<fcices  only  a  small,  almost 
insignificant  minimum  in  the  strain-depth  curve. 

(c)  There  is  a  secondary  strain  minimum  at  275  km. 

(d)  There  is  a  pronounced  minimum  at  400-475  km  depth  with  a  rapid 
increase  downwards  between  475  and  65O  km.  These  features  are  probably 
associated  with  the  depth  variation  of  pressure  and  temperature,  in 
combination  with  possible  phase  changes  or  changes  in  composition. 

(e)  The  strain  curve  is  almost  parallel  to  the  curve  showing  the 
variation  of  number  of  shocks  with  depth. 

(f)  The  average  strain  release  per  shock  is  practically  independent  0? 
depth,  the  decrease  f^und  being  statistically  insignificant. 
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